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\D ! ABSTRACT 

<N . 

Aims. This paper presents 452 new 21-cm neutral hydrogen line measurements carried out with the FORT receiver of the 
meridian transit Nangay radiotelescope (NRT) in the period April 2003 - March 2005. 

Methods. This observational programme is part of a larger project aiming at collecting an exhaustive and magnitude-complete 
[ HI extragalactic catalogue for Tully-Fisher applications (the so-called KLUN project, for Kinematics of the Local Universe 
studies, end in 2008). The whole on-line HI archive of the NRT contains today reduced Hi-profiles for 4500 spiral galaxies of 
' declination 5 > —40° I jhttp:/ /klun.obs-nancay.fr] |. 
Results. As an example of application, we use direct Tully-Fisher relation in three (JHK) bands in deriving distances to a large 
catalog of 3126 spiral galaxies distributed through the whole sky and sampling well the radial velocity range between and 
8 000 km s _1 . Thanks to an iterative method accounting for selection bias and smoothing effects, we show as a preliminary 
output a detailed and original map of the velocity field in the Local Universe. 
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This kind of HI data is crucial for constraining the gas 
and total mass function of spiral galaxies as a function of 
morphology and environnement, it allows also the map- 
ping of the total mass distribution from peculiar veloci- 
ties and thus provides strong constraints on cosmological 
models and large scale structure formation. They can in 
particular provide a unique starting point for total mass 
power spectrum studies. 

Study of peculiar velocities allows the verification of 
the current theory of cosmological structure formation by 
gravitational instability. It gives information on bulk mo- 
tion, and the value of f2 m (cf. reviews by Willick [20001 and 
Zaroubi 2002, and the comprehensive work by Strauss & 
Willick 1995). The velocity measurements are done using 
redshift independent secondary distance indicators, such 
as the Tully-Fisher (TF) relation for spiral galaxies, the 
Faber-Jackson, D n -a, Fundamental Plane (FP) relation, 
or the surface brightness fluctuations for early type galax- 
ies. The largest surveys so far are the Mathewson and 
Ford HHM1 sample, the MARK III (Willick et al. ITMTj) . 
SFI (Giovanelli et al. [EMI H aynes et al. 11955)1 . ENEAR 
(da Costa et al. I2000al2000b|) . and the updated FGC cat- 
alogue (2MFGC, Mitronova et al. l2004|) . Each contains in 
the order of one or two thousand independent distance es- 
timates in the local 80 h~ x Mpc volume. The use is then 
to compare them to the velocity field derived from the 
galaxy density distribution as infered from complete red- 
shift sample (e.g. PSCz, Saunders et al. 120001 or NOG, 
Marinoni et al. I1998[) 

Our own Kinematics of the Local Universe (KLUN) 
TF sample has been used in the study of Hq (Theureau 
et al. I1997) Ekholm et al. I1999|) and local structures 
(Hanski et al. ; 2001]) . The sample consists of all the galax- 
ies with published rotational velocities collected in the 
HYPERLEDAEI, database (Paturel et al. I2003ap . plus 
the recent large KLUN+ contribution (Theureau et al 
2005 and this paper). The total Tully-Fisher sample 
counts 15 700 spirals and uses five different wavelength 
galaxy magnitudes. B- and I-magnitudes come from vari- 
ous sources, carefully homogenized to a common system. 
The largest sources are DSS1 for B, and Mathewson et 
al. ([19921 [L996j) and DENIS (Paturel et al. [2004)) for I 
band. J, H, and K-magnitudes are from the 2MASS^| sur- 
vey (Jarret et al. |2"0U0|) . The 2MASS magnitudes, taken 
from a single survey, avoid any problems that the homog- 
enization may cause, and are thus exclusively used in data 
analysis. Further, we exclude the measurements with large 
errors and the galaxies that for other reasons, explained 
later in the text, are unsuitable for this study. 3126 galax- 
ies remain, which we use for the mapping of the peculiar 
velocity field within the radius of 80 Mpc. 

The paper is structured as follows: the Nangay radiote- 
lescope, the processing chain and the reduced HI data are 
presented in Sect. [2 The characteristics of the input Tully- 
Fisher catalogue are listed in Sect. [3l while the iterative 

5 http://leda.univ-lyonl.fr 

6 http:/ /www. ipac.caltech.edu/2mass/ 
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Fig. 1. Aitoff projection of the observed sample in J2000 
equatorial coordinates, histogram of radial velocities V20 
and HI fluxes (see Table [2]) 

method to obtain unbiased peculiar velocities from it is 
given in Sect. [4J finally, we give in Sect. [5] some prelimi- 
nary results and show some examples of peculiar velocity 
map realization. 

2. The HI data 

2.1. The Nangay observations 

The Nangay radiotelescope (France) is a single dish an- 
tenna with a collecting area of 6912 m 2 (200 x 34.56) 
equivalent to that of a 94 m-diameter parabolic dish. The 
half-power beam width at 21-cm is 3.6 arcmin (EW) x 22 
arcmin (NS) (at zero declination). The minimal system 
temperature at 6 = 15° is about 35 K in both horizontal 
and vertical polarizations. The spectrometer is a 8192- 
channel autocorrelator offering a maximal bandwidth of 
50 MHz. In this mode, and with two banks in vertical and 
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horizontal polarizations counting 4096 channels each, the 
spacing of the channels corresponds to 2.6 km.s -1 at 21 
cm. After boxcar smoothing the final resolution is typi- 
cally ~10 km.s" 1 . The 50 MHz bandwidth is centered on 
~ 1387 MHz, thus corresponding to an interval of 10,500 
km.s -1 centered on a velocity of 7000 km.s -1 (except for 
the few objects with a radial velocity known to be less than 
2000 km.s -1 , for which the observing band was centered 
on 5000 km.s -1 ). The relative gain of the antenna has 
been calibrated according to Fouque et al. (1990)); the fi- 
nal Hi-fluxes (Tabled are calibrated using as templates a 
set of well-defined radio continuum sources observed each 
month. 

One " observation" is a series of ON/OFF observational 
sequences; each sequence is made of ten elementary inte- 
grations of 4 seconds each, plus a set of 3 integrations 
of 2 seconds for the calibration, adding up in each cy- 
cle to 40+6 seconds for the source and 40+6 seconds for 
the comparison field. The comparison field is taken at ex- 
actly the same positions of the focal track as the source in 
the same cycle. In this way one minimizes efficiently the 
difference between ON and OFF total power. A typical 
meridian transit observation lasts about 35 minutes and 
is centered on the meridian, where the gain is known to 
be at its maximum; it contains a series of ^20 ON/OFF 
cycles. 

The processing chain consists of selecting good ele- 
mentary integrations or cycles, masking and interpolat- 
ing areas in the time-frequency plane, straightening the 
base-line by a polynomial fit (order in the range 1-6), and 
applying a boxcar smoothing. The maximum of the line 
profile is determined by eye as the mean value of the max- 
ima of its two horns after taking into account the rms 
noise (estimated in the base-line). The widths, measured 
at the standard levels of 20% and 50% of that maximum, 
correspond to the "distance" separating the two external 
points of the profile at these intensity levels. The signal 
to noise ratio is the maximum of the line (see above) over 
rms noise in the baseline fitted region. 

The total list of corrected Hl-astrophysical parameters 
(Tabic [2]), 21-cm line profiles (Fig. 3), and comments con- 
cerning the profiles (Table [3]) , are available in electronic 
form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
or via 



http: / / cdsweb.u-strasbg.fr / Abstract.html 



2.2. Sample characteristics and data description 

In the first five years of observations (2001-2005), since the 
upgrade of the Nangay receiver (FORT) , we have observed 
2500 galaxies, successfully detected about 1600 of them 
and fully reduced 1340 HI profiles. 

As a second KLUN+ release, we present here the spec- 
tra obtained for 452 of these galaxies, observed between 
April 2003 and March 2005. Some simple statistics is pre- 
sented on figure [2] The upper panel shows a comparison 
of some of our Hl-line width at the 20% level with equiv- 
alent measurements (Wp2o) found in the last large com- 



pilation of line widths by Springob et al. 120051 The over- 
lap is quite small, concerning 20 galaxies only. The few 
outlying galaxies marked with filled circles are identified 
either as distorted Hl-line, at the limit of detectability 
or Hi-confused with another neighbouring galaxy (case of 
pgc2350, pgc20363, pgc67934, pgc66850 and pgc54825, see 
Table |3|) . The other ones are well aligned on the first bi- 
secting line. Anyway, one could eventually guess a small 
systematic effect there : a slight over-estimation of the line 
width for large widths with respect to H Q or Arecibo mea- 
surement. This is explained easily by the general low signal 
to noise ratio we have for edge-on galaxies, in the range 
of fluxes we are concerned with. The Middle panel shows 
the distribution of signal to noise ratio S/N as a function 
of 20% level line width W20, and the bottom panel shows 
the rms noise a in mJy (outside the 21-cm line) versus 
integration time. In the latter, the curve shows the line 
a = 20/ VT~ t . 

Table [5] contains all the reduced HI parameters. Table 
[3] provides corresponding comments, when necessary, for 
each galaxy. Comments concern mainly object designa- 
tion, peculiar morphologies or peculiar HI line shape, spec- 
trum quality and HI confusions. The spectra and extracted 
data are assigned a quality code. A flag ' ?' or '*' warns for 
suspected or confirmed HI line confusion. The five quality 
classes are defined as follows: 

— A : high quality spectrum, high signal to noise and well 
defined HI profile 

— B : good signal to noise ratio, line border well defined, 
still suitable for Tully-Fisher applications 

— C : low signal to noise, noisy or asymmetrical profile, 
well detected but one should not trust the line width. 
The radial velocity is perfectly determined 

— D : low signal to noise, noisy profile at the limit of 
detection. Probably detected, but even radial velocity 
could be doubtful 

— E : not detected. The absence of detection, correspond- 
ing to the "E" code in the notes, is due to several pos- 
sible reasons: either the object was too faint in HI to 
be detected within a reasonable integration time (120 
full ON/OFF cycles, equivalent to 3 meridian tran- 
sit), which is probably the major case, or we did not 
know its radial velocity and it fell outside the frequency 
range, or the HI line was always behind a radar emis- 
sion or an interference... In a few cases, some standing 
waves are clearly visible in the full bandwidth plots (50 
MHz = 10,500 km.s -1 ). These are due to reflexions ei- 
ther in the cables, between the primary and secondary 
mirrors or between the secondary and tertiary mirrors. 
It happens when a strong radiosource (often the Sun) 
is close to the main beam of the antenna. Finally, when 
"no detection" is stated, the line was expected to fall 
within the observed frequency band and the value of 
the noise gives a fair upper limit for the HI signal. 

The distribution of the targets among the different 
classes is summarized in table [TJ 
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Table 1. Statistics of the detected galaxies vs. HI pro- 
file class. Among the "C" galaxies, 38 were flagged as HI 
confused 



profile class 


nb of galaxies 


A 


84 


B 


179 


C 


161 


D 


28 


E 


101 



The Aitoff projection of the catalogue in J2000 equa- 
torial is seen on Fig. f , together with the distribution of 
the radial velocities and HI fluxes. Most of the observed 
objects are in the range 4,000-10,000 km s^ 1 where the 
lack of Tully-Fisher measurements in the literature is the 
most critical. 



where Fjji is the observed raw Hi-flux, 
B f = ^(l + xT)(l + xt) 

T=(a 2 25 sin 2 /3 + b 2 5 cos 2 0)/6% w 
t=(a 2 25 cos 2 /? + fo 2 5 sin 2 (3)/6 2 NS 

0ew an d 9ns are the half-power beam dimensions of the 
Nangay antenna, (3 is the position angle of the galaxy de- 
fined north-eastwards, a 2 z and 625 are the photometric 
major and minor axis respectively. The parameter x is 
x = 0.72 ± 0.06 (Bottinelli et al.. [T990)) . The distribution 
of the East- West projection of D25 diameters is shown in 
Fig. 2. This is to be compared to the 4 arcmin width of 
the half-power beam. 



2.3. Data description 

Radial velocities Our observed radial velocities are listed 
in Table [2] (column 4) and correspond to the median point 
of the 21-cm line profile measured at 20% of maximum 
intensity. 

The internal mean error on V20 is calculated according 
to Fouque et al. (f990) as follows: 



S/N 



1/2 



where R is the actual spectral resolution, a = (W20 — 
W5o)/2 is the slope of the line profile, and S/N is the 
signal to noise ratio. The average of o"(Vao) is about 8 
km.s . 

Line widths Line widths are measured on the observed 
profile at two standard levels corresponding to 20% and 
50% of the maximum intensity of the line. The results 
listed in Table columns 6 and 9, have been corrected to 
the optical velocity scale. We also provide line widths cor- 
rected for resolution effect (Fouque et al ll990p in columns 
7 and 10. The mean measurement error is taken equal to 
3 ■ cr(V 20 ) and 2 ■ a(V 20 ) for the 20% and 50% widths, re- 
spectively. The data presented here are not corrected for 
internal velocity dispersion. Details about these correc- 
tions can be found in Bottinelli et al (1990), Fouque et al 
(PSD!) or in Paturel et al (|2003b|) . 



Hl-f luxes The detailed description of the flux calibration 
is given m Theureau et al. C2005). 

Hi-fluxes Fjfj (Table [H column 12) are expressed in 
Jy km s _1 . The values given in column 13 are corrected 
for beam-filling according to Paturel et al. {2003b ): 



F, 



HIc 



BfF, 



HI 



3. Building the Tully-Fisher sample 

3.1. The input data 

Rotational velocities Rotational velocities, i.e. the log V m pa- 
rameter used in the Tully-Fisher relation, have been mainly 
gathered from the HYPERLEDA compilation (16 666 galax- 
ies, Paturel et al. I2003a[) and were complemented by some of 
our own recent HI line measurements with the Nangay radiote- 
lescope (586 late type galaxies, Theureau et al. I2005| l. A few 
other measurements from Haynes et al.([1999J not previously in 
HYPERLEDA were also added. All our own 3300 HI spectra 
acquired with the Nangay radiotelescope antenna in the last 
decade were reviewed and assigned a quality code according 
to the shape of their 21-cm line profile. This study (Guilliard 
et al. 12004ft allows to flag efficiently several TF outliers due 
to morphological type mismatch or Hi-confusion in the ellon- 
gated beam of Nangay. Even if this Nangay subsample concerns 
only a part of the data (~ 25%), we substantially improved the 
apparent scatter of the TF relation. 

In the HYPERLEDA compilation, the log V m parameter is 
calculated from 2f-cm line widths at different levels and/or ro- 
tation curves (generally in H a ). The last compilation provides 
us with 50520 measurements of 2f-cm line widths or maximum 
rotation velocity. These data are characterized by some sec- 
ondary parameters: telescope, velocity resolution, level of the 
21-cm line width. For data homogenization, HYPERELEDA 
uses the so-called EPIDEMIC METHOD (Paturel, G. et al. 
2003bJ). One starts from a standard sample (a set of mea- 
surements giving a large and homogeneous sample : here, 
the Mathewson et al. 11996 data), all other measurements are 
grouped into homogeneous classes (for instance, the class of 
measurements made at a given level and obtained with a given 
resolution). The most populated class is cross-identified with 
the standard sample in order to establish the equation of con- 
version to the standard system. Then, the whole class is incor- 
porated into the standard sample. So, the standard sample is 
growing progressively. The conversion to the standard propa- 
gates like an epidemy. In summary, this kind of analysis con- 
sists in converting directly the widths for a given resolution r 
and given level 1 into a quantity which is homogeneous to twice 
the maximum rotation velocity (= log W, uncorrected for in- 
clination. A final correction is applied reference by reference to 
improve the homogenization. 
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Fig. 4. -excerpt. 21-cm line profiles of galaxies listed in Tabled profiles are classified according to their PGC name 
which is written above each panel. Ordinate and abscissae axes are graduated respectively in km s _1 and Jy. Note 
that heliocentric radial velocities are expressed in terms of optical redshift c-A^. The horizontal line represents the 
baseline of the profile, i.e. the zero flux level, from which the maximum is estimated. 
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The final log V m value is corrected for inclination : logV m = 
(log W — log (2 sin incl). Where the inclination incl is derived 
following RC3 (de Vaucouleurs [T99T|l : 

1 _ 1Q 2logR 25 

sin (incl) — — 

7?25 is the axis ratio in B at the isophot corresponding to 25 
magarcsec" 2 , logr = 0.43 + 0.053T for type T= 1 to 7 (Sa- 
Sd) and logr =0.38 for T=8 {Sdm). 

Magnitudes The 2MASS survey, carried out in the three in- 
frared bands J, H and K, collected photometric data for 1.65 
million galaxies with K s < 14 (Jarrett et al. [2000) and made 
the final extended source catalog recently public. Total mag- 
nitude uncertainties for the 2MASS extended objects are gen- 
erally better than 0.15 mag. We exclude any galaxy with the 
accuracy of magnitudes worse than 0.3. This accuracy appears 
reasonable when considering that it is almost impossible to ob- 
tain total magnitudes better than 0.1, due to the difficulty to 
extrapolate the profile in a reliable way. 

Extinction The extinction correction we applied includes a 
Galactic component, qg, adopted from Schlegel et al. (|1998[) . 
and a part due to the internal absorption of the observed 
galaxy, a\. Both depend on the wavelength. 

ttext = /g(A)cig + /i(A)ai. (1) 

The Galactic and internal wavelength conversion factors are 
/ G (A) =1.0, 0.45, 0.21, 0.13, 0.085 (Schlegel et al. [1998} and 
/i(A) = 1.0, 0.59, 0.47, 0.30, 0.15 (Tully et al. [T998l Watanabe 
et al. I2TT0T1 Masters et al. [213031 for B, I, J, H, and K bands, 
respectively. 

Radial velocities All radial velocities were taken from 
HYPERLEDA, often a mix and an average of several publi- 
cations and redshift surveys. Within the limit of 8 000 km s , 
we collected a sample of 32 545 galaxies. In general, when an HI 
measurement was available (i.e. for most of the TF subsample 
used in this study), the radio radial velocity was prefered, be- 
ing more acurate than the available optical one. As it is the use 
in HYPERLEDA, when more than two redshift measurements 
for the same galaxy were available, the most discrepant ones 
were rejected from the mean. The radial velocities are used at 
the first iteration of the IND method as a reference distance 
scale for the Malmquist bias correction. 

3.2. Selection and completeness 

The analytical treatment of the Malmquist bias effect with dis- 
tance, by applying the Iterative Normalized Distance method 
(IND), requires the strict completeness of the samples accord- 
ing to magnitude selection (Theureau et al. i 1998b] and Sect. 4.). 
The limits in magnitude are simply determined by eye as the 
knee observed in a logA^ vs. magnitude diagram, witnessing 
the departure from a homogeneous distribution in space with 
growing distance. This limit is in 'observed apparent magni- 
tude', independently of extinction or opacity correction. These 
corrections however are taken into account further, as part of 
the ND scheme itself, in what we call the effective magnitude 
limit (Sect. 4.1.). 



The adopted completeness limit are the followings : 
Jn m = 12.0, Hi im =11.5, and Ki im =11.0 (equivalent to Bi im ~ 
15 and In m ~ 13). Only the complete part of the sample in 
each band, about half of available data, is included for further 
study. 

The final selection is made according to the following con- 
ditions: 

— J, H, and K magnitudes < completeness limit 

— magnitude uncertainty < 0.3 mag 

— log of rotational velocity uncertainty < 0.03 

— T = 1-8 to keep only fair spiral galaxies 

— log i?25 > 0.07 to avoid face-on galaxies for which the ro- 
tational velocity is poorly determined 

After these restrictions there are 3263 spiral galaxies dis- 
tributed over the whole sky (see FigfS}. 

4. Method of analysis 

In this section we explain the Iterative Normalized Distance 
method for deriving the peculiar velocities. The 'iterative' 
means that a previously calculated peculiar velocity field is 
used for a more accurate estimation of new peculiar velocities. 
The 'normalized distance' is a quantity depending on the dis- 
tance and the absolute size or absolute magnitude of a galaxy, 
such that for any galaxy, the average selection bias (in the ter- 
minology of Strauss & Willick, I1995|) or the Malmquist bias of 
the second kind (according to Teerikorpi, il984[) can be given 
by a function depending on its normalized distance, the disper- 
sion of the distance criterion, and the completeness limit. This 
is illustrated in Fig. [6] where the TF residuals, plotted against 
the normalized distance modulus, clearly show the unbiased 
regime and the deviation due to the magnitude cutoff. 

A fully detailed description of the method follows, but we 
start by listing the main steps: 

1. Calculate the absolute magnitudes and the normalized dis- 
tances using the kinematical (redshift) distances. 

2. Calculate the TF relation using the unbiased part of the 
normalized distance diagram. 

3. Use the unbiased TF relations and the analytical 
Malmquist correction formula for estimating real space 
galaxy distances beyond the unbiased plateau limit (Fig5). 

4. Obtain the peculiar velocity field in a Cartesian grid in the 
redshift space by smoothing the individual peculiar veloc- 
ities given by the Malmquist corrected TF distances. 

5. Go back to step 1, and use the corrected kinematical dis- 
tances by substracting the smoothed peculiar velocity field 
values from the redshift velocities. 

This loop is repeated until converging values for the peculiar 
velocities are obtained. The peculiar velocities for all the galax- 
ies do not converge nicely, though. We thus extract the most 
unreliable galaxies (about 4 %) and recalculate the velocities 
with the reduced data set (see Sect. 4.4. and Fig. 6). As con- 
firmed by the tests done with a mock sample in Sect. 4.6. out- 
liers use to be mainly very low Galactitic latitude objects for 
which the corrected total magnitude is not well estimated from 
observed one, and core cluster members whose observed radial 
velocity, used as kinematical distance at the first iteration, does 
not reflect at all their true distance, due to their strong motion 
in the cluster potential. 
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Fig. 6. TF residuals (Y n against the normalized distance modulus (X n ). The "unbiased plateau" is the region at 
X n < — 1, where the binned average of the residuals, the balls with error bars, are close to the zero line. The curve 
going through the balls is the analytical correction formula described by Eq. 8 in Sect. 14.11 



4.1. Kinematical distances, normalized distances and 
absolute magnitudes 

Let us define the kinematical distance modulus as 

cz 

Mk = 51og — +25, (2) 
-no 

where cz is the observed heliocentric redshift, corrected by the 
CMB dipole motion. The absolute magnitude is 

M k = mo — (3) 

where m|j is the apparent magnitude, corrected for inclina- 
tion, extinction, and cosmological effects, as stated in Eq. 10 
of Paturel et al. (|1997[) (the cosmological correction is negligi- 
ble for all galaxies in the present study). 

If we consider that the TF relation is a linear law charac- 
terized by a given slope and a given dispersion (the zero-point 
being fixed either by some local calibrators or by adopting a 
value of Ho), and if we assume that the sample is actually 
complete up to a well defined apparent magnitude limit, then 
the selection bias at a fixed log V m is only a function of the 
distance (see Teerikorpi 1984] Theureau et al.[l997 1998b). In 
other words, the bias at a fixed log E m and at a given distance 
is only the consequence of the magnitude cut-off in the distri- 
bution of TF residuals and moreover it does not depend at all 
on space density law. 

By normalizing to a same luminosity class, i.e. a same 
log V m value, and by taking into account the variation of the ac- 
tual magnitude cut-off with extinction, one can build a unique 
diagram showing the bias evolution with distance. 

The distances and magnitudes are then scaled so that a 
sharp edge is seen at the sample completeness limit. 



The normalized distance modulus is defined as 

X n = /ik + (Mtf — mum + a oxt ) (4) 

where Mtf is the absolute magnitude of a galaxy, as given 
by the Tully-Fisher relation and mii m is the apparent magni- 
tude completeness limit. More explicitely, if we develop Mtf 
as [ai log Vm + ao] it appears that we normalize indeed to the 
same log Vm and the same effective magnitude limit mf^ = 
mum — a cxt . X n can also be expressed as Mtf — Mi im , i.e. 
the difference between the TF absolute magnitude at a given 
log Vm and the absolute magnitude cut off (in the TF residuals) 
at a given distance. 

The normalized magnitude 

Y n = M TF - A'/ k = Mtf -rn + ^ (5) 

corresponds to the departure of the absolute magnitudes cal- 
culated on the basis of kinematical distances from the true 
mean value given by the TF relation. This residual contains 
the contribution of magnitude and log E m measurement errors, 
internal TF dispersion, and peculiar velocity. 

Figure [S] shows normalized distance moduli vs. normalized 
magnitudes for the galaxies derived by the Tully-Fisher rela- 
tion. The curve going through the points is the analytical bias 
solution f(X n , cttf, Tliim) while the vertical line shows the up- 
per limit adopted for the unbiased normalized distance domain. 

4.2. Unbiased Tully-Fisher 

The TF relation states the linear connection between absolute 
magnitude and log E m 

M = ai log Vm + a (6) 
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One gets the slopes, ai, and zero-points, ao, of the TF re- 
lations in each band (IJHK) by least squares fit using the 
unbiased part of each sample. The unbiased part is the flat 
or plateau region in the normalized distance diagram. This 
subsample also provides us with the " 0- value" of the TF resid- 
uals that is used to compute the bias deviation for the whole 
magnitude complete sample (the horizontal line in FigHJ} . The 
method of deriving the TF relation in the unbiased plateau has 
been succesfully used in previous KLUN studies; a full statis- 
tical description can be found in Theureau et al. 11998 b and it 
was even numerically tested by Ekholm (|1997jl . 

In the iterative scheme one starts by assuming a priori val- 
ues for the TF slope and zero-point (here, only the slope is 
important and a rough value can be inferred directly from the 
whole sample). These values are used to compute the normal- 
ized distance and extract a first unbiased subsample. The loop 
"TF-slope — > normalized-distance — > unbiased-subsample — > 
TF-slope" can be repated a couple of times to be sure to start 
on the basis of unbiased values. 



4.3. Corrected distances 

The function f(X n , cjtf, miim) is deduced from the expectancy 
of Y n knowing X n assuming that the magnitude selection S(m) 
is described by the Heavyside function S(m) = 6(m — mn m ). 
We have then: 



E(Y n \X n ) = - 



where 



iii„ 



1 Y 2 
2-Ka 2cH 



Him = K(mi im ) = Mtf — m[ im + /ii, 

— MtF — (?Tllim — dext) + Mk 
- 1 2 



B = 



expl tt\cIY 

2-kg yl 2a 2 J 



From here we derive: 
f(Xtx,a,m\i m )=a\ — 



7r erfc(A) 



(8) 



where, 



A = 



and 



V2a 



X n 



erfc(x) — —= 



e y dy 



The corrected and unbiased distance modulus is then finally: 

fi c = ra - A/tf + f(X n , a, m Bm ) 

Note that Ho is cancelled out in Yu m : it is indeed hidden 
in the TF zero-point and explicitely present in fik but with an 
opposite sign. 

The reader will remark that our approach to the bias in this 
paper is radically different to what has been attempted with 
Marklll (Willick et al lT997l or Dekel et al [T999)) in which the ap- 
proach has been from the viewpoint of the classical Malmquist 
bias (using in particular some inhomogeneous density correc- 
tion). 



4.4. Iterations 

The peculiar velocities of galaxies are then smoothed onto a 
cartesian grid, v pcc (x,y, z) (Sect l4.5j) . As explained above, the 
method relies on kinematical distances; the normalized dis- 
tances and the absolute magnitudes needed for TF relation 
include the kinematical distance. These distances can be made 
more accurate by subtracting the peculiar velocities from the 
redshift: 



/i k = 5 log ■ 



Ho 



+ 25. 



(9) 



This distance modulus then replaces the value given by Eq. 
[2] and we can repeat the whole procedure with these updated 
distances. The new peculiar velocity field is again used for cor- 
recting the kinematical distances in the next iteration, and we 
keep on repeating the process until convergence. 

In practise, subtracting the whole v pec (x, y, z) as in Eq. 
[5] would overcorrect for the peculiar velocities, and cause di- 
verging oscillatory behaviour in the iterative process. Using a 
scaling factor A g (0.0, 1.0) so that 



5 log ■ 



>^v pec (x,y, 



Ho 



+ 25, 



(10) 



removes this problem. The superscript j corresponds to the it- 
eration number. Using the value A = 0.5 we reach converging 
values after about 5-10 iterations. Figure [7] shows this conver- 
gence for a few galaxies. Usually the v pec approaches nicely 
to a constant value, in some cases the values oscillate even 
with high j. We checked all these convergence curves by eye 
and rejected the worst cases. With the restricted sample we 
recalculated the peculiar velocities and used the results of the 
selected 3126 spirals galaxies for the peculiar velocity mapping. 



4.5. The tensor smoothing 

After deriving radial peculiar velocities of galaxies it is useful to 
interpolate these velocities at uniformly distributed grid points. 
The best method is to smooth the observed galaxy velocities 
with an appropriate window function. Dekel et al. (11999 [I dis- 
cuss the problems of smoothing a non-uniformly distributed 
set of radial velocities: 

The radial velocity vectors are not all pointing in the same 
direction over the smoothing window. Then, for example in a 
case of a pure spherical infall towards the window center, all the 
transverse velocities are observed as negative radial velocities. 
The net velocity in the smoothing window is then, incorrectly, 
negative, in stead of being zero. Dekel et al. {i!)!)9! call this 
the tensor window bias. They find that it can be reduced by 
introducing a local velocity field with extra parameters, which 
is to be fitted for the observed radial velocities in the smoothing 
window. The best results are obtained by constructing a three- 
dimensional velocity field with a shear, 



v(x) = B + L (x-Xc), 



(11) 



where L is a symmetric tensor, and x c is the window center. 
There are then nine free parameters, three for the actual win- 
dow center velocity B and the six components of the tensor L. 
The resulting grid point velocity is just v(x c ) = B. 

Furthermore, if the true velocity field has gradients within 
the effective smoothing window, a nonuniform sampling will 
cause an error, called the sampling-gradient bias. Dekel et al. 
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( 1999 ) suggest that this bias can be diminished by weighting 
the observed galaxy velocities by the volume V n , which is de- 
fined as the cube of the distance between the galaxy and its 
nth neighbour. This method gives more weight for galaxies in 
isolated areas. 

4.6. Testing the method 

We test these biases with a mock peculiar veloc- 
ity catalog. The mock catalog is constructed from 
the GIF consortium constrained n-body simulation of 
our 80/i _1 Mpc neighbourhood (Mathis et al. 120021 
http:/ /www. mpa-garching.mpg.de/NumCos/CR/ ). The 
simulation was run for a flat ACDM cosmological model, 
and it provides locations, velocities, masses, and luminosities, 
with and without the internal absorption, of 189 122 galaxies. 
The galaxy formation was defined by applying a semianalytic 
algorithm on the dark matter merger tree. We added the 
Galactic component of the absorption, as defined in Schlegel et 
al. ((1998 ), and selected the galaxies brighter than a magnitude 
limit. In the end there are 9800 galaxies with their apparent 
B band magnitude smaller than 14.5. 

Figure |S]shows the true vs. smoothed peculiar velocities us- 
ing different smoothing method and set of window parameters. 
These plots lead to the following comments: 

— the scatter is mainly related to the smoothing radius 

— outliers are essentially: 

- galaxies at low Galactic latitude, for which the magnitude 
is not well defined and where the sparser sampling leads the 
smoothing to diverge (vertical spreading) 

- galaxies belonging to clusters for which the kinematical 
distance derived from the redshift is strongly affected by 
the cluster velocity dispersion 

— one observes a tilt with respect to the line with slope = 
1, when no tensor window is used. This is the effect of 
the velocity gradient around structures and leads to an 
underestimation of infall patterns and to a cooler velocity 
field. 

The best \ 2 value is obtained for B + L smoothing, and 
Rsmooth = 0.15 x distance. This will be our choice of smoothing 
for the rest of the study. 

5. Preliminary results 

In this section we present the TF relation parameters obtained 
for three wavelengths (JHK) and some examples of maps of 
radial peculiar velocity fields, superimposed on the distribu- 
tion of galaxies. A more detailed kinematical study is beyond 
the scope of the present analysis and will be presented in a 
forthcoming paper. 

5.1. TF parameters 

Figures [§] show the final TF relations for the galaxies in the 
unbiased part of the normalized distance diagram. Tables [4]- 
list the parameters ai and ao and the scatter of the relation : 

M — ai log V m + ao (12) 

The observed scatter is comparable to what was found by 
Karachentsev et al 120021 using 2MASS magnitudes. The small 



Table 4. Tully-Fisher parameters: slope, zeropoint, scat- 
ter, and the number of unbiased plateau galaxies used in 
the relation. 





Oi 


a 


(Ttf 


N 


J 


-6.3 ± 0.31 


-8.8 ± 0.14 


0.46 


960 


H 


-6.4 ± 0.33 


-9.1 ± 0.14 


0.47 


1166 


K 


-6.6 ± 0.37 


-9.0 ± 0.16 


0.45 


990 



difference (we get slightly smaller cr's) can be explained easily 
by our optimization of the kinenatical distance scale through 
the iterative process described above. Accounting for the ob- 
served broadening due to apparent magnitude and log Vmax 
uncertainties and to the residual peculiar velocity dispersion 
affecting the kinematical distances, one obtains an internal 
scatter of ~ 0.44 mag for the TF relation in B and ~ 0.4 
mag in K. This is 0.1 mag greater than in studies restricted to 
pure rotation curve measurements of log Vmax ■ Here instead, a 
large majority of log Vmax measurement are from the width of 
global HI profiles. As we know, even once corrected for non- 
circular motions, this width is still determined by the shape of 
a galaxy's rotation curve, the distribution of HI gas in the disk 
and the possible presence of a warp fVerheiien ^OOlJ) . leading 
to a greater intrinsic Tully-Fisher scatter. 

5.2. Peculiar velocities 

Table [5] shows the first 10 KLUN+ galaxies with TF distances. 
The distances are expressed in km s _1 . That is followed by the 
kinematical distance, corrected by the radial peculiar velocity. 
This is our estimate of the true distance of a galaxy. Finally 
we list the observed redshift velocity in the CMB rest frame. 
The full catalog is given in electronic format only. 

Our peculiar velocities were obtained for all points in space 
having large enough galaxy density. We required that there 
should be more than 15 galaxies with peculiar velocity mea- 
surements within the smoothing radius around the point. Then 
we fit the 9-parameter tensor field to the peculiar velocities of 
these galaxies and set the value obtained at the center of the 
smoothing window (see Sect. 4 for more details). 

Since we use a distance dependent smoothing radius, the 
points close the Local Group must have a higher density of 
KLUN galaxies around them than the more distant points, 
for a succesful velocity field determination. This explains why 
some of the more distant grid points have a set value while 
there are apparently no galaxies around them. 

We compared the data to the Marklll distances (Willick et 
al. I1997p . Mark III catalog was compiled from six samples of 
TF and one of elliptical galaxies. It was converted to a common 
system by adjusting the zero points of the distance indicators. 
For the Malmquist bias correction the authors reconstructed 
the galaxy density field from the IRAS 1.2 Jy survey, and used 
it for the inhomogeneous correction formula. The corrected dis- 
tances for 2898 spirals and 544 ellipticals are publically avail- 
able and make a good comparison point for other peculiar ve- 
locity studies. 

Figure [TO] show the Malmquist corrected distances of in- 
dividual galaxies, as measured in Mark III, versus the corre- 
sponding value derived by us. The relative scatter dd/d is ~ 
0.2, correponding to an absolute uncertainty 0.43 in magni- 
tude scale, is fully compatible with the measured TF scatter 
(see Table|4]). A few points in Figure[l0]show a clear mismatch. 
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We found that these large discrepancies are due to errors in the 
input data in Mark III. These errors are listed and discussed 
in Appendix [X] 

Figure [11] shows the radial peculiar velocity field in the 
supergalactic plane, averaged over a disc having a thickness 
that increases towards the edge. The thickness of the disc in 
the center is zero, and its opening angle is 15°, so that at the 
edge (at 80 h _1 Mpc) the disk width is about 20 h _1 Mpc. 
The blue colors refer to regions where the radial peculiar ve- 
locity is towards us, the red regions are outfalling. The shade 
of the color corresponds to the amplitude of the motion, and 
is saturated at 1 000 km s _1 . The regions where the threshold 
requirement of at least two velocity measurements is not satis- 
fied are set white. The black dots are galaxies in HYPERLEDA 
with measured redshifts (not just the KLUN galaxies). Green 
circles mark some well known clusters. The maps are in real 
space coordinates, i.e. the redshift distances corrected by the 
peculiar velocity field. 

It is worth noting that in our v pcc maps one observes both 
the front and backside infall patterns around the main super- 
clusters and structures. It is particularly obvious on Fig. [11] 
for the regions of Virgo, Perseus-Pisces, N533, Norma, or even 
Coma, though it is located close to the limit of the sample. 
It seems that we even detect an outflow in the front side of 
the Great Wall. Similar features are seen on Fig |12l that shows 
other slices with different orientations in space. 

A wide region however, roughly centered on Centaurus 
cluster, seems to move away from us at a coherent speed of 
~400 kms" 1 on a scale greater than 20h~ 1 Mpc. The direc- 
tion and amplitude of this bulk motion are close to the one of 
the putative Great Attractor (Lynden-Bell et al 1988, Hudson 
et al 120041 Radburn-Smith 2006) and cannot be associated to 
any structure in particular. Anyway, it seems that this flow 
vanishes beyond a distance of 50/i _1 Mpc. 

As an example of quantitative result, we checked the am- 
plitude and direction of a bulk flow within a growing sphere 
centered on the Local Group. The result is shown on Fig. [13] 

At short scales the direction of the flow is compatible 
with most previous studies (Table [U Fig. 1131 bottom panel). 
In particular it coincides with the Great Attractor for R ~ 
20h~ 1 Mpc. At larger scales it first drifts towards the direction 
of the rich cluster region of Horologium-Reticulum, and after 
R ~ 40ft" 1 Mpc back to (I, b) ~ (310°, -8°). 

In the upper panel, we show the bulk flow amplitude within 
a sphere of growing radius (colored line). It oscillates strongly 
at short scales, as a consequence of the density heterogene- 
ity, and decreases to ~ 250 km s~ 1 at 40 h" 1 Mpc. Beyond 
this point it behaves more smoothly, as an indication that we 
reached the scale of the largest structures within our sample. It 
starts to rise (or oscillate) again beyond 60 h _1 Mpc, probably 
due to the sparser space sampling. 

The black solid curve shows the rms expected bulk velocity 
infered from the standard ACDM model: 



\rrms 



(v\R)) 



1/2 



2^2" 



P(k)W 2 (kR)dk 



1/2 



where P(k) is the mass fluctuation power spectrum and 
W (kR) is the Fourier transform of a top hat window of radius 
R. For the parametric form of P(k) in the linear regime, we 
use the general CDM model (see e.g. Silberman et al.. 12001) : 

p{k) = A c (n m , n A , n)T 2 {n m , a b , h- k)k n 



A c is the normalization factor and the transfer function T(k) 
is the one proposed by Sugiyama (1995). We restricted the 
analysis to a flat cosmological model with Q m +QA — 1, a scale- 
invariant power spectrum (n = 1), a baryonic density fib = 
0.04 (WMAP result, e.g. Spergel l20M]> . a Hubble constant fixed 
at h = 0.57 (which is our own preferred value, inferred directly 
from primary calibration by Theureau et al I1997|t1 . adjusting 
only the value of Q m . 

The best fit has been obtained for tt m = 0.22 in the dis- 
tance range 40-60 Mpc, where the value of {Vrms(R)) appears 
very smoothed. We also ploted the ± 0.02 curves around this 
best value. What we observe confirms WMAP results on fl m 
(e.g. Spergel [2006J and seems coherent with the expected rms 
bulk velocity within a sphere for standard ACDM model (see 
e.g. Willick [20001 or Zaroubi [2002]) . thus with no bulk motion. 

One should be prudent anyway in such kind of conclusion 
: the theoretical prediction is here the rms value of a quantity 
that exhibits a Maxwell distribution (see e.g. Strauss 199TJ); a 
single measurement of the flow field is only one realization out 
of this distribution and gives only very weak constraints on the 
cosmological model. 
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Appendix A: Mark III errata & rejections 

When closely inspecting the Mark III data we found a few 
inaccuracies. In comparison we used the Mark III catalog pro- 
cat. 



vided by the CDS archives, http://cdsweb.u-strasbg.fr 



VH/198, and the data given by the HYPERLEDA database, 
http://leda.univ-lyonl.fr, as they were presented in May 2003. 
A few values in Mark III were replaced by those listed in 
HYPERLEDA. Some of the Mark III galaxies were rejected, 
having large differences to the HYPERLEDA values. 

Table [O lists galaxies having their PGC numbers incor- 
rectly identified in Mark III. We list first the number given in 
Mark III, followed by the correct number, alternative name, 
and the Mark III data set including the galaxy. 

Table lA~2l lists galaxies that were rejected due to their sus- 
picious values for redshift velocities (values given in the CMB 
rest frame). All the galaxies with \czm3 — czhl\ > 200 km s _1 
were studied, for most we maintained the Mark III values. 

Table IA.3I has the galaxies rejected due to the log V m 
uncertainties. Here we considered galaxies with | log Kn,M3 — 
logVn^HLl > 0.15. Notice that Mark III logW^ values are here 
converted to the inclination corrected log Vm of HYPERLEDA. 
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Fig. 3. Projection of the major axis D25 on the East- West 
direction 

KLUN + 




Fig. 5. Sky distribution of KLUN galaxies used in the 
current analysis 
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Fig. 2. Top : comparison of some of our HI- line width 
at the 20% level with some independent measurements 
from Springob et al. [2Q05 : . Middle : distribution of signal 
to noise ratio S/N as a function of 20% level line width 
W2o- Bottom : rms noise a in mJy (outside the 21-cm 
line) versus integration time. The curve shows the line 
g = 20/ VT~ t 
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Fig. 7. The evolution of peculiar velocities with respect 
to the number of iterative loop. Here are the two first 
galaxies in our sample, and an exceptionally "bad" case 
(pgc 17978) which was rejected from the final analysis. 
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Fig. 9. Tully-Fisher relations in J, H, and K bands for 
unbiased plateau galaxies. See Table 2] for the parameters 
of the relations. 
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Fig. 10. Comparison of Marklll and KLUN distances ex- 
pressed in km.s -1 . The colors represent the different Mark 
III samples (Willick et al. [T39T|) : HMCL is black, W91CL 
red, W91PP green, CF blue, MAT turquoise, and A82 
purple. 
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Fig. 13. Upper panel shows the amplitude of the observed 
bulk motion within a growing sphere of radius i? ma x, and 
the corresponding ACDM expected curve for a flat cos- 
mological model using : n = 1, fi& = 0.04, h — 0.57, and 
ft m = 0.22 ± 0.02 . The bottom panel show the direction 
of the bulk flow, in galactic coordinates, and its evolution 
with i? max = lO-60/i" 1 Mpc. The main published results 
as listed in Table [6] are also shown. 
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pgc/leda 


A T 4 TV f n 

NAME 


RA (2000) DEC 


V 2 o 


a v 


W20 \V20c 


0" 14/20 


W50 


W 50c 


<y wso 


F(HI) F(HI) C 




S/N 


rms 


Q flag 


pgc0000115 


PGC000115 


J000145. 3-042049 


























2.2 


E 




pgc0000287 


NGC7813 


J000409. 1-115902 


9122. 


13. 


403. 


396. 


39. 


387. 


385. 


26. 


1.84 


1.85 


0.51 


3.1 


1.4 


B 




pgc0000317 


ESO409-011 


J000439. 3-282738 


8041. 


4. 


291. 


284. 


11. 


282. 


280. 


8. 


4.08 


4.12 


0.68 


7.9 


2.8 


A 




pgc0000432 


ESO409-016 


J000556. 1-310610 


7837. 


7. 


335. 


328. 


21. 


320. 


318. 


14. 


3.46 


3.48 


0.74 


5.6 


2.8 


A 




pgcOOOlOll 


NGC0054 


J001507.7-070625 


5333. 


18. 


429. 


422. 


55. 


413. 


411. 


37. 


3.17 


3.23 


1.22 


2.2 


3.2 


B 




pgc0001165 


PGC001165 


J001759.4-091620 


6979. 


22. 


440. 


433. 


67. 


418. 


416. 


45. 


2.49 


2.54 


0.98 


2.1 


2.5 


C 




pgc0001431 


PGC001431 


J002221.4-012045 


4898. 


6. 


173. 


166. 


17. 


77. 


75. 


11. 


7.47 


7.52 


0.75 


17.2 


4.1 


C 


c 


pgc0001453 


PGC001453 


J002238. 2-240733 


9951. 


18. 


544. 


537. 


55. 


504. 


502. 


36. 


3.09 


3.09 


0.75 


3.5 


2.0 


c 


c 


pgc0001542 


NGC0102 


J002436. 5-135722 


7333. 


26. 


464. 


457. 


77. 


438. 


436. 


51. 


1.64 


1.66 


0.67 


2.0 


1.7 


c 




pgc0001813 


NGC0131 


J002938. 5-331535 


1422. 


25. 


205. 


198. 


75. 


183. 


181. 


50. 


10.19 


10.39 


6.87 


1.9 


27.3 


c 




pgc0001897 


PGC001897 


J003058. 2-092415 


























1.5 


E 




pgc0002046 


ESO540-002 


J003414. 1-212812 


6983. 


5. 


287. 


280. 


14. 


275. 


273. 


9. 


3.17 


3.20 


0.51 


7.4 


1.9 


A 




pgc0002047 


ESO540-001 


J003413. 6-212619 


8048. 


3. 


131. 


124. 


9. 


92. 


90. 


6. 


4.21 


4.33 


0.35 


22.0 


2.0 


A 




pgc0002164 


PGC002164 


J003611. 3-323428 


4381. 


16. 


202. 


195. 


48. 


162. 


160. 


32. 


1.93 


1.95 


0.63 


4.0 


2.6 


C 


c 


pgc0002333 


PGC002333 


J003920. 1-102855 


10948. 


11. 


328. 


321. 


32. 


319. 


317. 


22. 


1.82 


1.84 


0.69 


2.8 


2.3 


B 




pgc0002349 


NGC0178 


J003908. 5-141018 


1447. 


3. 


138. 


131. 


8. 


99. 


97. 


5. 


9.82 


9.99 


0.74 


22.9 


4.0 


A 




pgc0002352 


NGC0192 


J003913.4+005151 


4212. 


11. 


423. 


416. 


32. 


390. 


388. 


21. 


2.27 


2.46 


0.43 


5.5 


1.3 


C 


c 


pgc0002369 


ESO540-009 


J003921.0-185503 


3893. 


10. 


188. 


181. 


29. 


172. 


170. 


19. 


1.39 


1.42 


0.46 


4.2 


2.0 


B 




pgc0002380 


NGC0187 


J003930. 2-143922 


3937. 


9. 


319. 


312. 


28. 


293. 


291. 


19. 


4.29 


4.33 


0.81 


5.5 


2.6 


B 




pgc0002424 


PGC002424 


J004029. 1-101820 


8104. 


10. 


336. 


329. 


31. 


317. 


315. 


21. 


1.49 


1.50 


0.38 


4.2 


1.3 


B 




pgc0002460 


UGC00435 


J004059. 6-013802 


5438. 


5. 


299. 


292. 


15. 


281. 


279. 


10. 


2.57 


2.62 


0.47 


8.8 


2.3 


C 


c 


pgc0002476 


PGC002476 


J004122.4-354836 


6408. 


9. 


258. 


251. 


28. 


240. 


238. 


19. 


2.10 


2.12 


0.60 


4.6 


2.5 


B 




pgc0002637 


IC1578 


J004426. 0-250434 


6737. 


34. 


390. 


383. 


102. 


306. 


304. 


68. 


1.79 


1.80 


0.58 


2.7 


1.6 


C 




pgc0002797 


ESO411-016 


J004739.6-275651 


1783. 


9. 


156. 


149. 


27. 


120. 


118. 


18. 


2.38 


2.39 


0.57 


6.7 


2.9 


A 




pgcuuuzyoo 


iioU4tt-Uzz 


J 005042. (-312302 


5959. 


1 n 
19. 


145. 


138. 


06. 


no 

y2. 


nn 
90. 


OO 

38. 


3.47 


3.54 


1.1 1 


n 
3.9 


inn 
10.9 


C 


? 


pgc0003636 


UGC00627 


J010100.6+132806 


11757. 




461. 


454. 




468. 


466. 




2.06 


2.08 


0.66 


2.6 


1.7 


C 




pgc0003642 


ESO351-031 


J010060.0-351433 


























2.3 


E 




pgc0004316 


ESO352-018 


J011209.4-321432 


9893. 


11. 


246. 


239. 


34. 


237. 


235. 


22. 


1.43 


1.47 


0.64 


2.7 


2.4 


C 


c 


pgc0004352 


PGC004352 


J011235. 1-040824 


5664. 


32. 


341. 


334. 


97. 


318. 


316. 


64. 


1.11 


1.13 


0.71 


1.5 


2.1 


D 




pgc0004641 


ESO352-034 


J011725. 2-354702 


9628. 


11. 


277. 


270. 


34. 


228. 


226. 


23. 


2.68 


2.72 


0.55 


6.2 


2.2 


B 




pgc0004703 


ESO542-003 


J011844.5-193736 


6444. 


23. 


423. 


416. 


68. 


374. 


372. 


45. 


1.34 


1.37 


0.46 


3.1 


1.5 


C 




pgc0005055 


UGC00928 


J012312.4-003828 


























2.3 


E 





Table 2. Astrophysical Hi-parameters (excerpt). 

Column 1: PGC or LEDA galaxy name; 

Column 2: most usual galaxy name; 

Column 3: J2000 equatorial coordinates; 

Column 4: systemic heliocentric radial velocity (km s _1 ); 

Column 5: rms error (km s -1 ); 

Column 6: total line width at 20% of the maximum intensity (km s _1 ); 
Column 7: total corrected line width at 20% (km s _1 ); 
Column 8: rms error (km s" 1 ); 

Column 9: total line width at 50% of the maximum intensity (km s _1 ); 

Column 10: total corrected line width at 50% (km s" 1 ); 

Column 11: rms error (km s _1 ); 

Column 12: observed Hl-fiux (Jy km s _1 ); 

Column 13: beam-filling corrected Hl-flux (Jy km s _1 ); 

Column 14: rms error (Jy km s _1 ); 

Column 15: signal to noise ratio; 

Column 16: rms noise; 

Column 17: quality code (see Sect. 2) 

Column 18: flag ("c" indicates confirmed HI confusion with the emission of another galaxy; "?" means that confusion 
is suspected but not certain) 
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pgc/leda 


Type 


log D 2 5 


P.A. 


Q 


comments 






log(O.l') 


deg. 






pgc0000115 


Sab 


0.87 


92.5 


E 


small galaxy group, also pgc 131 (MCG-01-01-023) close to the [...] 


pgc0000287 


Sb 


0.91 


158.0 


B 


Theureau et al 2005 


pgc0000317 


SO-a 


0.96 


37.0 


A 




pgc0000432 


Sa 


1.03 


8.7 


A 




pgcOOOlOll 


SBa 


1.14 


92.6 


B 


Paturel et al 2003, pgcl024961 also in beam 3' NW, late type, [...] 


pgc0001165 


SO-a 


0.96 


126.3 


C 


SO galaxy 


pgc0001431 


Sbc 


0.74 


48.0 


Cc 


HI confusion w UGC212 at V=4840, multiple/interacting galaxy, [...] 


pgc0001453 


111 


0.35 




Cc 


HI confusion w ESO473-018 at V=9923 Theureau et al 2005 


pgc0001542 


SO-a 


0.97 


126.5 


C 


SO-a B galaxy 


pgc0001813 


SBb 


1.25 


62.5 


c 


=NGC131, edge of NGC134 m cf at V=1579 but prob no HI confusion [...J 


pgc0001897 


Sb 


0.91 


123.5 


E 




pgc0002046 


SBb 


0.98 


125.4 


A 


=ESO540-002, galaxy group, pgc2047 and pgc2057 in the beam at [...] 


pgc0002047 


SBc 


1.12 


168.6 


A 


pgc2057 (early type) also in the beam at V=8172 


pgc0002164 


SO-a 


0.95 


40.8 


Cc 


HI confusion w ESO350-037=pgc2157 at V=4312 


pgc0002333 


Sab 


0.91 


174.0 


B 




pgc0002349 


SBm 


1.33 


3.0 


A 


Bottinelli et al 1982, interacting or peculiar galaxy 


pgc0002352 


SBa 


1.31 


164.9 


Cc 


HI confusion, galaxy group in the beam w NGC197,NGC196, [...] 


pgc0002369 


Sb 


1.04 


172.2 


B 




pgc0002380 


SBc 


1.17 


149.2 


B 




pgc0002424 


Sab 


0.85 


65.0 


B 




pgc0002460 


Sab 


0.98 


24.8 


Cc 


probably HI confusion w pgc090496 in cf, V unknown 


pgc0002476 


SBbc 


0.95 


56.0 


B 




pgc0002637 


Sb 


0.93 


17.7 


c 


Theureau et al 2005 


pgc0002797 


Sab 


1.01 


46.5 


A 




pgc0002955 


SBb 


1.10 


168.5 


C? 


ESO411-021 Irr in beam 3' N, unknown V, interaction ? probably [...] 


pgc0003636 


Sa 


0.88 


118.5 


c 




pgc0003642 


SBb 


0.93 


174.3 


E 


successfully observed by Theureau et al 2005 


pgc0004316 


Sa 


1.04 


128.1 


Cc 


confusion with NGC0427 (Sa, V=10012) at the edge of the beam 11' [...] 


pgc0004352 


Sa 


1.05 


89.0 


D 


lenticular ? low SNR, pgc4346 E galaxy V=5611 also in beam 


pgcUUU4641 


SBb 


0.98 


171.6 


13 


multiple object ? HI spectrum OK 


pgc0004703 


SO-a 


1.00 


62.4 


c 


low SNR 


pgc0005055 


SO 


1.01 


47.7 


E 


t=-1.8 


pgcOOG-5112 


SO-a 


0.98 


164.4 


E 


t=-0.5 barred galaxy 


pgc0005145 


SO-a 


1.02 


64.0 


E 


t=0.0 


pgc0005151 


Sb 


1.05 


30.8 


D 


Theureau et al 2005 


pgc0005168 


SBb 


0.76 


107.7 


C 


Theureau et al 2005 


pgc0005505 


Sa 


0.94 


72.0 


B 


pgc737697 V=9232 at edge of beam 3' SE 



Table 3. Notes on Hi-observations (excerpt). 

Column 1: PGC or LEDA galaxy name 

Column 2: morphological type from HYPERLEDA 

Column 3: logarithm of isophotal D25 diameter in 0.1 arcmin from HYPERLEDA 

Column 4: Major axis position angle (North Eastwards) from HYPERLEDA 

Column 5: quality code and Hi-confusion flag "c" (confirmed) or "?" (possible) (see Sect. [3]) 

Column 6: comments; conf="HI confusion", comp=" companion" , cf=" comparison field", poss=" possible" ,w="witlr' 
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Fig. 8. Shown here are the true peculiar velocities (on y axis) vs. the velocities obtained with a smoothing method for 
the GIF consortium simulated data. The upper left corner of each panel gives the essential information about the plot; 
TS stands for the tensorial smoothing, using nine parameters, the bulk velocity B and the tensor L, BK is the three 
parameter bulk smoothing (only B), and SC means scalar smoothing, using a simple Gaussian smoothing window. 
r s is the smoothing radius, either a distance dependent value, e.g. 0.15<i or a fixed value, e.g. 10 Mpc hr 1 . The rms 
difference between the smoothed and the true velocity field, 6v rms , is given in km s _1 . N is the number of points, 
the simulated galaxies, in each figure. Upper left figure shows the smoothing giving the smallest Sv ims . In each figure 
the galaxies in densest regions were excluded, as well as the objects at the low Galactic latitudes. The two panels 
in the upper right corner show how these points give substantially more divergent values. The "Dense regions" and 
< 15°" maps were obtained with the TS r s — 0.15d smoothing. V4 and Vs correspond to the volume weighted 
method, using the volume defined by the fourth and the eight closest galaxy, respectively. The volume weighting is 
claimed to reduce the sampling gradient bias, but in this test they fail to produce a better estimate of the true velocity 
field than the non-weighted method. 
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Table 5. An excerpt of the table showing the distance data. We list the name and the galactic coordinates of the 
galaxies, followed by the three TF distances. Then we give the corrected kinematical distance, g?£ = — V pcc . All 
are expressed in km s _1 . The last entry is the redshift velocity, corrected to the CMB rest frame. 





I 


b 


^TF.J 






d c k 


U 3k 


PGC 2 


113.96 


-14.70 


6084 


6030 


5992 


4762 


4751 


PGC 54 


109.57 


-33.17 


8864 


8936 


9152 


8771 


8379 


PGC 76 


109.81 


-32.67 


7997 


7592 




6968 


6570 


PGC 102 


111.34 


-27.22 


5214 


5125 


5211 


4834 


4720 


PGC 112 


110.61 


-30.23 


4961 


4754 


4764 


4493 


4449 


PGC 120 


108.40 


-37.98 


3893 


3771 


3815 


4202 


4051 


PGC 129 


108.42 


-37.97 


6413 


6155 


6105 


4177 


4026 


PGC 176 


94.33 


-63.84 


7662 


7359 




6412 


6117 


PGC 186 


107.24 


-42.50 


9047 


8557 


9158 


8034 


7541 


PGC 195 


355.68 


-77.39 


6281 


5862 


6319 


6254 


6538 



Rot around X 9= 0.0° 





Fig. 11. Radial peculiar velocity field in the supergalactic plane. Blue colors are regions with negative peculiar veloc- 
ities, red colors refer to positive ones. Black dots are all galaxies with a known redshift. Green marks show positions 
of some well known structures - here Virgo, Ursa Major, Norma, Pisces, N533, and Zwicky33. The coordinates are in 
"real space", i.e. redshift distances corrected for the smoothed peculiar velocity field, in units of km s _1 . 
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Fig. 12. The peculiar velocity maps, as in Fig. [TTJ projected on four discs of different orientations with respect to 
the supergalactic plane (SGZ=0). Figure a is rotated by 10°around the SGX axis. It illustrates well the features 
around two large rosary structures: the Great Wall starting from Coma and extending counterclockwise up to the 
Centaurus region, and another starting from Perseus-Pisces complex, going down through N533 and beyond. Both 
front and backside infall are visible all along the structures. Figure b is rotated by -30°around SGZ from the SGY=0 
plane, and shows the large infall motions towards the two opposite regions, one in Perseus-Pisces and the other in 
the Great Attractor or Centaurus area. Figure c, rotated 50°around SGY from the SGX=0 plane, shows the infall 
patterns towards Pegasus, Hydra, and Coma. Figure d, that is almost perpendicular to the supergalactic plane, 
rotated -30°around SGY from the SGX=0 plane, shows an example of some very detailed structure of the velocity 
field. 
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Table 6. Directions, in Galactic coordinates, of some of the main bulk flow measurements or large galaxy concentra- 
tions. The number corresponding to Figll3l is given in parentheses. 





I 


b 


ref. 


SNIa (1) 


282 


-8 


Riess et al. 119971 


ENEAR (2) 


304 


25 


da Costa et al. I2000al 


SFI+SCI+II (3) 


295 


25 


Dale & Giovanelli 120001 


SMAC (4) 


260 


-1 


Hudson et al. 119991 


PSCz (5) 


260 


30 


Rowan-Robinson et al. 2000 


LP 


343 


52 


Lauer & Postman 119941 


CMB 


276 


30 


Local GrouD motion, Koeut et al. 119931 


GA 


307 


9 


Great Attractor, Lvnden-Bell et al. 119881 


HR 


270 


-55 


Horoloaium-Reticulum. Lucev et al. 119831 


SC 


315 


29 


Shapley Concentration, Scaramella et al. [1989 



Table A.l. Mark III errata: PGC numbers. 



Mark III 


correct 


alt. name 


Mark III data sets 


PGC 10631 1 


PGC 95735 




HMCL, W91CL, CF 


PGC 64575 


PGC 64632 


NGC 6902 


HMCL 


PGC 57053 


PGC 57058 


UGC 10186 


W91CL 


PGC 71291 


PGC 71292 


UGC 12572 


W91CL 



1) PGC 10631 = UGC 2285 is projected on PGC 95735, but at lower redshift. Mark III clearly refers to the latter, but lists it 
as UGC 2285. 

Table A. 2. Mark III errata: velocities. 



Mark III HYPERLEDA Mark III data sets 



PGC 26680 


7266 


12736 


HMCL 


PGC 72301 


11747 


7348 


W91PP 


PGC 17136 1 


4417 


8870 


CF 


PGC 26561 2 


2136 


1664 


MAT 


PGC 67258 


2278 


2684 


MAT 


PGC 62889 


2513 


5745 


MAT 


PGC 30753 3 


2946 


3813 


MAT 


PGC 9551 


4057 


4633 


MAT 


PGC 20324 4 


4114 


5681 


MAT 


PGC 15790 


4202 


6199 


MAT 


PGC 47832 


4477 


4827 


MAT 


PGC 31723 5 


4866 


4130 


MAT 


PGC 62411 6 


5075 


5996 


MAT 


PGC 5964 


5859 


5406 


MAT 


PGC 3144 


5974 


5053 


MAT 


PGC 8888 


6070 


5049 


MAT 


PGC 64523 


6384 


5271 


MAT 


PGC 19363 


7092 


6666 


MAT 


PGC 2001 


7112 


6180 


MAT 


PGC 44349 


9696 


9916 


MAT 



1) Value incorrectly copied from the data source 

2) Lower value suggested by three independent sources 

3) See note at Giovanelli et al. [19971 (Antlia 27146) 

4) Possible confusion with galaxy 1' SE 

5) A typographical error in MAT file? See Mathewson et al. 119921 Fig. 3 

6) Unclear Ha observation 
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Table A. 3. Mark III errata: logV B 





Mark III 


HYPERLEDA 


Mark III data sets 


PGC 72024 


2.321 


2.131 


HMCL,W91CL,W91PP,CF 


PGC 47707 1 ' 2 


2.139 


2.391 


HMCL 


PGC 65338 1 


2.079 


1.844 


HMCL 


PGC 7111 1 


1.978 


2.251 


HMCL 


PGC 36875 1 


2.253 


2.453 


W91CL 


PGC 9841 3 


2.049 


1.596 


W91PP 


PGC 51784 


2.209 


1.889 


CF 


PGC 5453 


2.429 


1.980 


CF 


PGC 17113 1 


1.994 


1.830 


MAT 


PGC 16201 


1.807 


2.092 


MAT 


PGC 32821 1 


2.023 


1.861 


MAT 


PGC 67818 1 


2.052 


1.843 


MAT 


PGC 39139 


1.760 


2.018 


MAT 


PGC 51982 1 


2.284 


2.057 


MAT 


PGC 24328 


2.232 


2.013 


MAT 


PGC 64042 1 


2.442 


2.200 


MAT 


PGC 13778 1 


2.295 


2.123 


MAT 


PGC 2383 


1.978 


2.150 


MAT 


PGC 36875 1 


2.241 


2.453 
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1) These galaxies have large differences in inclination stated in Mark III and in HYPERLEDA. 

2) Mark III states the inclination as 90°, while it is 34° in HYPERLEDA. The latter is probably correct. 

3) Several sources favour the HYPERLEDA value. 



